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Abstract 

A method of building and investigation of the Fermi surfaces for three-dimensional 
crystals subjected to a uniform magnetic field is presented. The Hamiltonian of a 
charged particle in the crystal is treated in the framework of the zero-range poten- 
tial theory. The dispersion relation for the Hamiltonian is obtained in an explicit 
form. 
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1 Introduction 

The shape of the Fermi surface determines the kinetic and equilibrium properties of the 
electron gas in the crystal matter £Q as well as the dynamics of a single electron in crystal 
The shape of this surface can be reconstructed drastically by a high magnetic field 
j3IUE]. However, a uniform magnetic field changes totally the translation properties of an 
electron in the crystal lattice: the appearance of a new length scale (the magnetic length) 
leads to the famous phenomena related to the "commensurability-incommensurability" 
transitions. In particular, a fractal structure arises in the diagram describing the depen- 
dence of the two-dimansional electron spectrum on the magnetic flux (Azbel'-Hofstadter 
butterfly jHlED- It was shown recently that this fractal structure is an inherent charac- 
teristic not only of two-dimensional spectral diagrams but also of the "flux-energy" and 
"angle-energy" diagrams for a three-dimensional magneto-Bloch electron 0E]- 

The translation symmetry of the Bloch electron in a uniform magnetic field is de- 
fined by the magnetic translation group ^U|, which has more complicated structure in 
comparison with the translation group without the field. In particular, the ordinary 
quasi-momentum is no longer a conserved quantity, therefore, at high magnetic fields a 
modification of the definition of the Fermi surface is required. Such a modification is given 
by V. Ya. Demikhovskii with coworkers fUE]- ^ n P a P ers [TT1 112j . the construction and 
investigation of the Fermi surfaces in the magnetic Brillouin zone have been performed 
by means of the tight-binding approximation. In particular, it was shown that the change 
of the topology of the Fermi surface under the influence of the magnetic field can cause 
the metal-semiconductor transition. 



The procedure of obtaining the tight-binding Schrodinger operator (also called the 
Harper operator) leads to a rather rough approximation of the initial Hamiltonian; for 
example, in the 2D case it is known that the "flux-energy" diagram for the original 
periodic Landau operator differs essentially from the classical Hofstadter butterfly [T3*l ll4|. 

In this connection, it is interesting to build and investigate the Fermi surfaces in the 
magnetic Brillouin zone for the periodic Landau operator without additional approxima- 
tions. Here we consider the case of the Landau operator perturbed by a lattice A of 
short-range scatterers. To get an explicitly solvable model we deal here with a limiting 
case of the short-range potential, namely, with the zero-range potential E3 El > also 
called "point potential". We stress, that the method of zero-range perturbations for ob- 
taining the Fermi surfaces in 3D crystals without magnetic fields was successfully used 
earlier |18| . On the other hand, the spectral properties of the three-dimensional Landau 
operator perturbed by a periodic point potential was studied in papers [HJ HHJ |20|. In 
particular, in the paper [20J, a version of the magneto-Bloch analysis was proposed which 
makes it possible to get an explicit form of the dispersion relation for the problem we 
consider here. 



2 The model Hamiltonian and dispersion relation 

We start with the quantum-mechanical Hamiltonian H of a single electron subject to a 
uniform magnetic field B 

ff ° = ^(p->>) 2 > w 



2/x v c 

where e and n are the charge and mass of the electron, respectively, and A(r) is the vector 
potential of the field B. We choose the symmetric gauge: A(r) = B x r/2; in this case 
the form of Hq in any Cartesian coordinate system with the z-axis directed along the field 
B is independent of the choice of the other coordinate axes. 

We use the standard notations uj = \eB\/cLi (the cyclotron frequency), Im = (H/liuj) 1 ^ 2 
(the magnetic length), and O = 2irhc/\e\ (the magnetic flux quantum). Throughout the 
paper we denote by b, b = — T$/<f>o, the vector of the magnetic flux density: if ai and a 2 
are nonzero vectors, then |b(ai x a2)| is the total number of the magnetic flux quanta 
through the parallelogram spanned by ai and a 2 . Note that |b| = {2'k1 2 m )~ 1 . 

All important physical information concerning the Hamiltonian H is contained in the 
Green function Go at energy E defined as Go(r,r'; E) = (r\(E — Hq)^ 1 ^'). It is known 
(see, e.g., Ref. (2D 122]) that this function can be represented in the following forms: 

G (r, r'; E) = $(r, r^r - r'; E) = $(r, r')F 2 (r - r'; E) , (2) 

where 
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Here rn is the projection of r on the direction of the field B and r_|_ = r — r||; the function 
Li(x) is the l-th Laguerre polynomial. 

It is convenient to decompose G into a sum of singular and regular (with respect to 
the limit r — > r') parts : 

G,(r, r'; E) = ^ W X ^ + Gg«(,, r'; E) . (6) 

Note that G r e9 (r, r; E) is well defined at r = r' and 

were ({s,v) is the generalized Riemann zeta-function (in other words, the Hurwitz zeta- 
function) [23] . 

Now we consider a crystalline lattice T with nodes 7; the Bravais lattice of the crystal 
we denote by A. Choose a basis ai, a 2 , a 3 of A, then a basic cell Cr of T (the Wigner- 
Seitz cell) is fixed: Cr = {tiai + £2^2 + t3a 3 : < U < 1}. If K is the set of nodes of Y 
containing in Cr, then T = K + A; this means that each vector 7 G T can be represented 
in the form 7 = k + A, where k G K, A G A (note that such a representation is unique). 
Without loss of generality we will suppose that G K. We choose the potential V(r) of 
the crystalline lattice V in the form 

KG K AgA 

where V re (r) is the confinement potential near the node k,. To obtain an explicitly solvable 
model, we pass to the zero-range limit for V K . This means that we consider V K as the 
limit of potentials of the form c K W(r — n), where W(r) ~ outside a small sphere of 
radius R centered at zero, the coupling constant c K is of order R, and J W{r) dr = 1. At 
the limit R — > the potential V n is characterized by one parameter only, namely, by the 
scattering length p K , which is related to the binding energy of the ground state for V K by 

E K = 

(see Ref. [IH])- What is important, at the zero-range limit the Green function G(r, r'; E) 
of H has the following explicit expression in terms of the Green function Go{r, r';E) of 
H (see Ref. US EH): 

G(r,r';£) = G (r, r';£)- £ G (r, 7; E) (S^^E))^^, G (-f r'; E), (9) 

7,7' er 

where S~ l (E) is the matrix inverse to the infinite matrix S(E) with elements 

5^, + (l-5^,)Goh,~/';E). (10) 



The Hamiltonian H = Ho + V in some sense is a three-dimensional analogue of the one- 
dimensional Kronig-Penney Hamiltonian. The spectral properties of H can be extracted 



from the properties of invertibility of the infinite matrix S(E). Namely, at least for 
E < hu/2 the energy E belongs to the spectrum of H if the matrix S(E) is not invertible. 
Under so-called "rationality condition" ^U|j we can reduce the problem to invert the 
infinite matrix S(E) to a problem of finite-dimensional algebra. According to JU], we 
call the field B rational with respect to A if for some basis ai, a 2 , a 3 of A, all the 
numbers b(aj x a^) are rational. It is clear that every B obeys this condition after 
an appropriate infinitely small change of direction, therefore, the rationality condition 
imposes no essential restrictions on B. If B is rational with respect to A, then a basis slj 
of A can be chosen in such a way that r] = b(ai x a2) > and b(a2 x a 3 ) = b(ai x a 3 ) = 0. 
Denote rj = N/M, where N and M are coprime positive integers; in our approach the 
crucial role is played by the following (MK) x (MK) matrix S, where K is the number 
of nodes in K: 

Sq >q '(p,E) = exp[— mm'bfa' x a 2 )] x 

oo 

^ 5(Aiai + (A 2 M + m)a 2 + k, m'a 2 + k', A 3 a 3 ; E) x 

Ai,A2,A3=— oo 

exp [7ri(Aia! + (A 2 M + m)a 2 )(b x k - r)\ 1 (M\ 2 + m) - 2A • p)] . (11) 

Here q denotes the pair (m, n) with k & K, m — 0, . . . ,M — 1 and p is the magnetic 
quasimomentum. The magnetic Brillouin zone is defined by the inequalities < p\ < 
M -1 , < p 2 < 1, < p 3 < 1 and, in this zone, the dispersion relation for H reads 

det S(p, E) = . (12) 

For a fixed p, equation (JT2)l has infinitely many solutions E s (p), s = 0, 1, . . .; each E s (p) 
is an M-fold degenerate eigenvalue of H. Dispersion laws E = E s (p) are continuous 
functions of p; the set of all values of E s (p) is an interval J s . In a natural way, the 
intervals J s form clusters called magnetic bands or Landau bands, each magnetic band 
contain KM adjacent intervals J s (subbands) . For the Fermi energy Ep, the equation 

E s (p) = E F (13) 

determines a periodic surface in the quasimomentum space (the Fermi surface). 

3 Results of numerical analysis 

Using Eqs. (fTT|) — (JT3j) we analyze numerically the shape of the Fermi surfaces in a crystal 
with the simple-cubic lattice for various directions of the field B. The lattice constant 
a = 7.5 nm is chosen relevant to the geometric parameters for the 3D regimented quantum 
dot superlattice considered recently in Ref. [21]. As to the scattering length, we put 
p ~ 1 nm, this corresponds to the binding energy E ~ 30 meV. 

Let the magnetic field B is directed along the edge a 3 of the cubic elementary cell and 
7] = 1; this value of 77 corresponds to the field strength ~ 70 T (note in this connection 
that it was reported on the possibility to achieve the magnetic field strength 28 MG 
|25j). In the considered case, the lowest magnetic band of the spectrum has the form 
[— 1.21£q; — 0.68eo], where e — fw/2. The Fermi surfaces for various energy values from 



this band are depicted in Fig. 1. The warped ellipsoids in the panels (a) and (d) are the 
electron and hole surfaces respectively; in this case there exist only closed trajectories of 
the charged particle in the p-space. Other topologies of the Fermi surface are shown in 
panels (b) and (c); according to the definition from [2] these surfaces have the topology 
of range 2 and 3 respectively. 

Surprisingly an inclination of B to the edge of the cubic cell causes a profound recon- 
struction of the Fermi surfaces. As an example we consider B = 5(0,3/5,4/5). In the 
new basis &[ = a 1; a' 2 = a 2 + a 3 , a' 3 = 3a 2 + 4a 3 , we have a 3 ||B and rj = 1/5. Therefore, 
the magnetic Brilluoin zone is defined by the inequalities < p\ < 1/5, < P2,P3 < 1, 
and the lowest magnetic band splits into 5 subbands. Fig. 2 shows variation of the Fermi 
surface as the Fermi energy Ep varies in the first subband [— 1.11 £o! — 0.98 eo]- The elec- 
tron and hole Fermi surfaces have the shapes of warped cylinders (panels (a) and (d) 
respectively) or highly twisted sheets (panels (b) and (c)). It is particularly remarkable 
that in this case the surfaces are open in the pi-direction. In going from the first sub- 
band of the lowest magnetic band to the higher subbands, the shape of the Fermi surface 
changes; examples of this change are shown on Fig. 3, where the Fermi energy Ep lies in 
the second and fourth subbands respectively. 

In the examples above, the filed B is parallel to a face of the cubic cell. We have built 
the Fermi surface for the field of the form B = 5(12/25, 9/25, 4/25), where r] = 1/25 but 
by the technical reasons do not give the corresponding figure here. In this case, the Fermi 
surface is divided into great numbers of connected pieces having the shape of winding 
sheets. Such a complex structure of the Fermi surface reflects the fractal structure of the 
spectrum of the periodic Landau operator . 
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